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Abstract 

Photometric observations of GO Cyg were performed during the July-October 2002, 
in B and V bands of Johnson system. Based on Wilson's model, the light curve analysis 
were carried out to find the photometric elements of the system. The 0-C diagram which 
is based on new observed times of minima suggests a negative rate of period variation 
(4r < ^) ^'^^ system. 

1 Introduction 

GO Cygni (HD196628; a : 20'*37'^,5 : +35°26') has been classified as a /3-Lyr type eclipsing 
binary of short period of about P = 0'^.71 whose visual magnitude varies from a maximum 
light value of Vmax = +8.6 to a minimum light value of Vmin = +9.25 and has a period 
variation. The results of the light curve analysis and the period study can be found in many 
papers such as Ovenden (1954), Popper (1957), Rovithis et al.(1990), Hall & Louth (1990), 
Sezer et al.(1993), Edalati & Atighi (1997), Rovithis et al.(1997), and Oh k Ra (1998). 

The structure of the paper is as follows. In Sec. 2 we introduce the process of observations 
and find the times of minima. In Sec. 3 the period study of the system is presented. The 
Wilson-Devinney model has been used to get the orbital and physical parameters of the 
system which are given in Sec. 4. In Sec. 5, the absolute dimensions of the system is calculated 
by the use of spectroscopic data together with the parameters estimated from the present 
study. A brief discussion and our main conclusions are presented in Sec. 6. 



2 Observations and times of minima 

We present the observed light-curve of GO Cygni, together with the new times of minima from 
the photometric observations carried out at Biruni Observatory. The 20 inches cassegrainian 
reflector which is equipped with an uncooled RCA4509 multiplier phototube, was used for 
differential photometry in B(A = 4200A) and V(A = 5500A) filters of Johnson's system. The 
observations were performed during Aug. /Sep. 2002 versus HD 197292 and HD 197346 as 
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the comparison and the check stars. In each observing run in B or V filter, the integration 
time was fixed for 10 seconds. The observed data were transferred to computer using an A/D 
converter. The reduction of data and atmospheric corrections were done using REDWIP 
code developed by G. P. McCook to obtain the complete light-curves in the two filters (Fig. 
1). The times of minima in Heliocenteric Julian Date (two primary and two secondary) were 
calculated by fitting a Lorentzian function to the observed minima data points (Dariush et 
al. 2003). Fig. 2, shows a sample Lorentzian fit to the observed primary minimum of GO 
Gyg on Aug. 31 in filter B. 

3 Period study 

Though a constant period for GO Cyg was first measured by Szczyrbak (1932) other studies 
such as Purgathofer & Prochazka (1967), Cester et al. (1979), Sczer et al. (1985), Jones et 
al. (1994), Rovithis et al.(1997), and Edalati & Atighi (1997) suggested a period variation 
in GO Cyg. For all the observed photoelectric times of primary and secondary minima from 
1951 to 2002 the computed (O-C)'s versus epoch is plotted in Fig. 3. For each minima, the 
residuals are calculated according to the ephemeris given by Hall &; Louth (1990) 

Min.I = JD2433930.40561 + 0*^.71776382 x E. (1) 
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Figure 2: A sample Lorentzian fit to the observed primary minimum data on Aug. 31 in 
filter B. 



In equation 1, £^ is the number of cycles. Since the period is changing with time an extra 
quadratic term can be added to equation 1 and rewrite it as 



Min.I = J£)2433930.40561 + 0''.71776382 xE + ^xE"^ 



(2) 



Prom the method of least square the coefficient of quadratic term in equation 2 is easily found 
to be ^ = 9.35 x 10~^^. As it is seen from Fig. 3a, the 0-C residuals suggest a polynomial 
function in the form of 

n 

ATiE) = J2cjE^, (3) 

where Ar(£^) is the differences between the observed and calculated times of minima for any 
cycle E. In the case of GO Cyg, the 3^^ order polynomials {j = 3 in equation 3) seems to 
be fitted better than the parabola (j = 2 in equation 3) to the observed times of minima. 
The value of the polynomial coefficients cj and the rms error of the approximation are listed 
in Table 1. The residuals of the observed 0-C from AT{E) which do not deviate more than 
ibO.Old are shown in Fig. 3b and 3c. Using equation 3, the real period of the GO Cyg and 
its relative rate of change can be calculate according to the following formula (see Kalimeris 
et al. 1994) 

P{E) =Pe + AT{E) - AT{E - 1), (4) 

n-l 



dP 



^^^^ = dE = + + - 1)'}' 



(5) 



j=0 



j=0 
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Coefficient 



Second 



Third 



Co 

Cl 
C2 
C3 

rms error 



-2.92771 X 10-^ 
-6.38377 X 10-'^ 
1.26577 X 10-^° 

0.00625 



1.45 X 10-^ 
-2.3711 X 10-^ 
3.01298 X 10"^° 
-4.41579 X 10-^5 
0.00597 



Tabic 1: Coefficients of second and third order approximation of the (0-C) diagram of GO 
Cyg. 



where Pg is an ephemeris period given in equation 1. Fig. 4 shows the period of GO Cyg as a 

function of time and its relative rate of change according to the 3'''^ order polynomial fitting 
function on the 0-C diagram. From Fig. 4, it is obvious that at present time, the period is 
decreasing. We refer the reader to Sec. 6 where we presented a completed discussion about 
the period variations. 

4 Light curve analysis 

The photometric solution and light curve analysis is done using the Wilson and Devinney's 
(1971) model which is based on Roche model. Before doing analysis we binned the observed 
data in B and V filters presented in Fig. 1 in equal phase intervals and applied Wilson's LC 
code in mode 4. This mode is for semidetached binaries with star 1 accurately filling its 
limiting lobe, which is the classical Roche lobe for synchronous rotation and a circular orbit, 
but is different from the Roche lobe for nonsynchronous and eccentric cases. The applied 
constraints arc that fii has the lobe filling value and that L2 is coupled to the temperatures. 
The LC code consists of a main FORTRAN program for generating light curves. Solutions 
were obtained for each of the B and V filters, separately. The preliminary values for the 
orbital and physical elements of the system were adopted from the previous studies. The 
third light and eccentricity were supposed to be equal to zero (e = 0, Z3 = 0). In order to 
optimize the parameters, an auxiliary computer programme was written to compute the sum 
of squares of the residuals (SSR) in both B and V filters, using LC outputs. The set of 
parameters {i,q,T2,^i,^2, Li, and L2) given in Table 2, are derived after optimization in 
the way as we described. Fig. 5 and 6 show the optimized theoretical light curves together 
with the observed light curves in B and V filters. These theoretical curves, correspond to the 
optimized parameters given in Table 2. In Table 4 we tabulated our results (as an average in 
B and V filters) together with previous photometric solutions from other sources. 

5 Absolute dimensions 

By combining the photometric and spectroscopic data, it is possible to estimate the absolute 
physical and orbital parameters of the system. The first spectroscopic elements of the system 
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Figure 3: (a) Plot of (O-C)'s versus epoch for all the observed photoelectric times of primary 
and secondary minima from 1951 to 2002. The residuals are calculated according to equation 
1. The dotted line and solid line represent a 3'"'^ order polynomial (j=3 in equation 3) and a 
parabola (j=2 in equation 3) fitting functions to all of the residuals respectively (see Table 
1). (b) Residuals between the observed 0-C differences and the best fitted parabola, (c) 
Residuals between the observed 0-C differences and the best fitted 3^*^ polynomial. 
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Figure 4: The period of GO Cyg as a function of time (continues line) and its relative rate 
of change (dotted line) according to the equations 4 and 5 with j=3. The difference P — Pg 
is referred to the ephemeris period Pe=0-71776382d which is given in equation 1. 



6 



Parameter 

i (degree) 

Ti (fixed) 

T2 

Ai 

A2 

91 

92 

x\ 

Li 
L1+L2 

L2 
L1+L2 

ri (pole) 
ri (point) 

ri (side) 
ri (back) 
r2 (pole) 
r2 (point) 
r2 (side) 
r2 (back) 
S(0 - C)2 



Filter B 

4200 

77.72±0.02 

0.428±0.001 

2.734 

2.905±0.005 

10350 
6528±30 



Filter V 

5500 

78.38±0.02 

0.428ib0.001 

2.741 

2.902±0.002 

10350 
6404±40 



1 


1 


0.5 


0.5 


1 


1 


0.32 


0.32 


0.82 


0.58 


0.79 


0.51 


0.965 


0.947 


0.035 


0.053 


0.427 


0.427 


0.586 


0.586 


0.455 


0.455 


0.482 


0.482 


0.259 


0.259 


0.300 


0.300 


0.267 


0.267 


0.285 


0.285 


0.00172 


0.00214 



Table 2: Optimized parameters of GO Cygni. 
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Figure 5: Theoretical light curve (LC output; solid line) in filter B according to the param- 
eters given in Table 2. Points are the observational data. 



were prepared by Pearce (1933) and a more recent spectroscopic orbit was obtained in 1988 by 
Holmgren using the cross-correlation techniques ( Sezer et al., 1992). We use the spectroscopic 
elements of Holmgren together with the orbital parameters and physical parameters which 
have been calculated already from our observation to calculate the absolute dimensions of GO 
Cyg. The results of these calculations based on equations 6 (see Kallrath & Milone, 1999) is 
given in Table 3. In equations 6, /(Mi, M2, i) is mass function, Mtot = Mi + M2 is the total 
mass, Ki is the radial velocity amplitude of component 1, a is the relative semi-major axis, 
and g is the acceleration due to gravity at the surface. 

^ f{Mi,M2,i)(l+q? 

/(Mi,M2,i 



M2 _ „( Ml 
Mq ~ yvM0 

_ M/Mq . 
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Figure 6: Same as Fig. 5 but in V filter 

6 Results and discussion 

A closer look at Fig. 5 and 6, shows that the primary minimum is symmetric around phase 1 
and the secondary minimum occurs at phase 0.5. Also, the out of eclipse portion of the light 
curve in two figures is the same which is consistent with the results of Sezer et al. (1985) and 
Jassur h Puladi (1995). But there is a lack of symmetry in the secondary minimum of both 
figures around phase interval of 0.6 to 0.7 which was also reported earlier (Sezer et al., 1993 
and Edalati &: Atighi, 1997). Since the 0-C diagram of GO Cyg clearly shows a variation in 
period ,the asymmetry in secondary minimum is most probably due to gas streaming or mass 
transfer. The result of our study support this idea that the primary component M\ has filled 
its Roche lobe, and the second component (M2; less massive and cooler one) is somewhat 



Parameter 


Filter B 


Filter V 


A(A) 


4200 


5500 


Ml (Mo) 


2.74 


2.82 


M2(M0) 


1.73 


1.75 


Ri{Rq) 


2.16 


2.17 


R2iRQ) 


1.09 


1.08 




58.72 


57.54 


L2{Lq) 


2.17 


2.13 



Table 3: Absolute elements of GO Cygni. 
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smaller than its Roche lobe. Let the total mass of the binary system be Mt^t = Mi + M2 
and take the total angular momentum of the orbiting pair as J orb- In reality, the total 
angular momentum Jtot is composed of two terms J orb and Jspin(spin or rotational angular 
momentum) such that Jtot = Jorb + Jspin but in practice the contribution of the total spin 
angular momentum is no more than 1-2% of Jorfe, so it can be ignored (Hilditch, 2001). Now 
consider a conservative mass transfer between the two components as well as the conservation 
of total orbital angular momentum Jgrb- Therefore M = ^^^^ = or dM\ = —dM2 (since 
all the mass lost by one component is gained by its companion), also ^^3^ = 0. Using the 
relation between Mtot and Jgrb together with Kepler's third law, lead us to the following 
relation (Hilditch, 2001) 

3P(Mi - M2) ^ ^ 

in which P and P are -in order- the orbital period and the rate at which the orbital period 

is changing and Mi is mass of the massive (and hotter) component. 

Now suppose that the system follows path 1 on the 0-C diagram which is defined by the 
parabolic fit. In this case according to the equations 4 and 5, P > (also Mi > M2), so 
it is clear from equation 7 that Mi > and the mass of the more massive component (Mi) 
must increases with time. But from the Roche model wc expect Mi < since the primary 
component has filled its Roche lobe and is loosing its mass to its cooler companion M2. As 
a result if we consider a simple conservative model for mass transfer with the spin angular 
momentum ignored then there is an inconsistency between the rate of period change and 
conservative model of mass transfer of the system. This deviation may be due to the role of 
stellar winds or spin-orbital angular momentum interchange. 

On the other hand if the system follows path 2 on the 0-C diagram where is defined 
by a 'S^'^ order polynomial, then it is clear from Fig. 4 that at the present time and after 
1996, P < 0, which means that equation 7 is valid and the mass of the primary component 
is transferring to the secondary one {Mi < and M2 > 0) which is in agreement with the 
Roche model we considered for the GO Cyg. The point is that the ascending branches of 
both parabolic and S''^ order polynomial fit to the 0-C data points in Fig. 3a intersect 
around period cycle £^=24000 (see also Rovithis et al. 1997), so the observations after this 
time is very important to give a clear understanding of the period change of GO Cyg and its 
behaviour. The results of our observations strongly support this idea that the rate of period 
change is negative but we must extend our observational baseline in order to collect more 
minima to study the behaviour of the GO Cyg 0-C diagram and to sec which model is most 
probable. If the new times of minima still follow path 2 in 0-C diagram then there is no 
doubt that the increase in rate of change of the period of the system has been stopped and 
the system is now undergoing a decrease of period. 
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Parameters 


Ovenden 


Mannino 


Popper 


Sezer et al. 


Rovithis et al. 


Edalati et al. 


Present 




(1954) 


(1963) 


(1980) 


(1993) 


(1997)1 


(1997) 


study^ 


z (degree) 


78.777 


79.880 


77.606 


79.722 


78.414 


82.311 


78.05 


Ti 


10350 


10350 


10350 


10350 


10350 


10350 


10350 


T2 


5605 


5904 


5966 


6043 


6912 


6428 


6452 


Li 
L1+L2 


0.969 


0.961 


0.972 


0.954 


0.931 


0.931 


0.956 


L2 
L1+L2 


0.031 


0.039 


0.028 


0.034 


0.069 


0.068 


0.044 


ri (back) 


0.482 


0.482 


0.482 


0.482 


0.555 


0.481 


0.482 


ri (side) 


0.455 


0.455 


0.455 


0.455 


0.504 


0.455 


0.455 


ri (pole) 


0.427 


0.427 


0.427 


0.427 


0.445 


0.426 


0.427 


r2 (back) 


0.282 


0.271 


0.221 


0.278 


0.303 


0.275 


0.285 


r2 (side) 


0.265 


0.256 


0.214 


0.261 


0.280 


0.258 


0.267 


r2(pole) 


0.257 


0.249 


0.210 


0.254 


0.265 


0.251 


0.259 



Table 4: Photometric solutions from other light curves. 

1. Based on the observations in Bucharest. 

2. Average on B and V filters. 
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